A PRIORI ESTIMATES OF SMOOTHNESS OF 
SOLUTIONS TO DIFFERENCE BELLMAN 
EQUATIONS WITH LINEAR AND QUASILINEAR 

OPERATORS 

N.V. KRYLOV 

Abstract. A priori estimates for finite-difference approximations 
for the first and second order derivatives are obtained for solutions 
of parabolic equations described in the title. 



1. Introduction 

The goal of this article is to prove a priori estimates for solutions of 
finite-difference approximations of parabolic Bellman equations with 
linear and quasilinear operators. In the latter case the nonlinear oper- 
ator defining the equation is still supposed to be convex with respect to 
the second-order derivatives of the unknown function. We present esti- 
mates for the finite-difference approximations of the first and second or- 
der spatial derivatives. In particular, our results cover finite-difference 
approximations for degenerate quasilinear parabolic equations. As far 
as we are aware these are the first results for such equations. The main 
parts of the linear and quasilinear operators entering Bellman equa- 
tions are assumed to be linear a^Afc operators, that is written as a lin- 
ear combination of pure second order derivatives in certain directions 
that are common to all operators. This assumption is always satisfied 
if the equation is uniformly nondegenerate and is generally necessary 
if we want to restrict ourselves to monotone difference approximations 
and meshes that are obtained from a fixed one by scaling (see more 
about it in Remark 12.41 below) . Our results are valid for usual Bellman 
equations and also for optimal stopping and impulse control problems 
associated with them. 

The motivation to obtain a priori estimates is the following. There 
is an approach suggested in [TU], [H], and to establishing the rate 
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of convergence of Uh to u as h \, 0, where u is the true solution, Uh 
the solution of finite-difference approximation of the same equation, 
and h typically is the mesh size. Two main ideas of this approach 
are that the original equation and its finite-difference approximation 
should play symmetric roles and that one can "shake the coefficients" 
of the equation in order to be able to mollify under the sign of nonlinear 
operator. 

For elliptic Bellman equations with constant coefficients and Lips- 
chitz free terms the first idea led to the rate of convergence of order 
^1/3^ for generic finite-difference approximations and h^^"^ in the case of 
OfcAfc operators (see Remark 1.4 and Theorem 5.1 in [TD], also see [2]). 
In contrast with the popular belief that assuming more smoothness of 
the data does not lead to better rates of convergence, it is proved in [5] 
that if the free terms are in C^'^, then the rate is at least h for constant 
coefficient a^Afc case and h"^ for equations with better structure. 

The second idea was introduced to treat equations with variable 
coefficients and led to quite satisfactory error bounds for u — Uh from an 
"easy" side (depending on how the equation is written this can be either 
upper or lower estimate of u — Uh). To get an estimate from the other 
side on the basis of the idea of symmetry between the approximating 
and the original equations one needed to solve the following problem: 

(P) in the case of variable coefficients estimate how much the solu- 
tion of the finite-difference equation loses in the process of shaking the 
equation. 

In the absence of solution of the problem (P) the idea of symmetry 
was still useful but only in obtaining some intermediate estimates (see, 
for instance, [2] and pLlj) and various approaches to getting the error 
bounds from the "hard" side were developed. In addition to the above 
cited papers the interested reader should consult [3], [1], and the ref- 
erences therein. Note that for generic finite-difference approximations, 
under the assumptions of Theorem 5.3 of [llj the result of [4J is the 
same /i^/^, but the result of Theorem 5.4 of [TT] is improved from /i^/^^ 
to h^'\ The issue of solving the problem (P) for generic finite-difference 
approximation remains unsettled and it is not clear how far off /i^^^ is 
from the true rate. 

The problem (P) was recently reduced to the problem of estimating 
the modulus of continuity of approximate solutions and solved in [13] 
for OfcAfc case in which a sharp error bound of order h^/"^ was obtained. 
The idea of symmetry worked again as in the constant coefficients case. 
This activity was continued in [7], where for the first time equations 
in domains were treated, and in [6j, where under various smoothness 
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assumption the rates /i^/^, h, and h'^ are obtained for linear degener- 
ate equations of ctfcAfc form. For linear case the rate h^^"^ was earher 
obtained in [8j by a method close to a method from [12] (Lemma 5.1 
of [8| is a version of Theorem 2.1 of [12] )■ However, this method does 
not allow one to get rates h and h!^. 

The main technical result of [13] is the a priori estimate of the de- 
rivative of Uh with respect to x stated as Theorem 5.2 and proved by 
quite subtle estimates. It turns out that there is a much easier method 
to prove Theorem 5.2 of [13] which in addition carries over to much 
more general equations with quasilinear operators and to obtaining es- 
timates for the second-order finite differences of Uh- The method is 
almost as simple as the one used in [6] for linear equations. 

We present this new method here and concentrate only on a priori 
estimates to keep the article within reasonable limits. Once the a priori 
estimates are obtained, one can follow familiar patterns to get error 
bounds in various cases of linear or quasilinear operators, degenerate 
or weakly nondegenerate or else uniformly nondegenerate, with or 
C^'^ coefficients. In particular, we hope to obtain first estimates on 
the rate of convergence in the case of degenerate quasilinear operators 
with Lipschitz coefficients. Our preliminary computations also show 
that under the assumptions of Theorem 12.121 the estimate \u — Uh\ < 
jY^2/3 j-^Q^g the elliptic case. These and some other indicated below 
possible applications of our results we intend to develop in the future. 

Hongjie Dong and the referees of the paper made valuable comments 
on the first version of it for which the author is sincerely grateful. 

2. Setting and main results 
Our first few results concern equations of the type 

F{6'^u, Ah/^u, Sh/^u, u) =0, (2. 1) 

where 

^(0, qk, Pk,i^) = F{(t), q, p, ijj, t, x) 

di 

= sup[r"(t)0 + V (a^(t, x)qk + bt{t, x)pk) - e{t, x)ij + T (p, V^, t, x)], 
|fc|=i 

(2.2) 

Ah/^u, Sh,ef,u are finite-difference approximations of the time de- 
rivative, the pure second-order derivative in direction £k, and the first- 
order derivative in direction ik, respectively. Detailed description of 
the above objects now follow. 
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Let A be a separable metric space, d,di > 1 integers, and let 

r" = r-(t), at = at{t,x), 6- = 6^(t,x), c" = c^{t,x), 

be real-valued bounded functions of {a,t,x) defined on A x M x M*^ 
for k = ±1, ...,±di. Also let some vectors ik G be defined for 
k = ±1, ±di and let 

r,/io e (0,oo), 5g(0,1], Ko,KuK2,K3 e[0,oo), meR (2.3) 

be some constants fixed throughout the article. It is worth noting that 
ik, k = ±1, ...,±di, are not supposed to form a basis in M'^ or even 
generate W^. This becomes crucial when one proves the estimates of 
the first-order differences of solutions with respect to parameters on 
which the coefficients may depend. Notice also that the lengths of ikS 
can be different and some of them can be just zero (and we will use 
this possibility later). The constant T gives us the time interval [0,T), 
on which the equation is investigated, ho "calibrates" the mesh-sizes 
in X variable, the constant 6 will appear in various requirements of 
nondegeneracy. The constant Kq is the most basic one, it is used in 
formulations of the very basic assumptions. The constant Ki is used to 
control either the maximum magnitude of the solution or its oscillation. 
The constant K2 will appear in our assumption on the growth of / with 
respect to the "gradient" of the solution (see Assumption 12. 51 (ii). which 
looks very much like the one commonly used in the theory of quasilinear 
PDEs. By the way, the author's efforts to use Assumption 12.51 (iii), 
stated similarly, failed.) The constant is used to control various 
quantities having lesser impact on our results than those controlled by 
Kq, Ki, K2. Finally, the constant m is used to extract various results, 
which in the theory of parabolic PDEs one gets after replacing u{t, x) 
with u{t, x)e"^K 

For any vector / G M'^, r/, r > 0, and function u introduce 

/ N u(x + ril)—u(x) ,, 

V 

. . u(t + TT(t),x) - u(t,x) . , , u(t + T,x) - u(t,x) 

d^U[t,X) = , drU{t,X) = , 

r r 

, , u{x + Tjl) — 2u{x) + u{x — Tjl) 

where the notation = (l/2)(|a| ± a) is used. Observe that with the 
above definition of 6j equation (12. ip makes perfect sense for t < T for 
functions u{t,x) defined only for t <T. We do not need to extend u 
beyond T in order to compute the finite-difference approximation of its 
derivative in time for t < T. 
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Assumption 2.1. (i) The functions r", a^, 6^, and c° are continuous 
with respect to a; 

(ii) the functions 6^ satisfy the Lipschitz condition with constant Kq 
with respect to x; 

(iii) the function satisfies the Lipschitz condition with constant 
with respect to x; 

(iv) we have 

i-k = -ik, a% = al \ik\<Ko, r°>0, > 
(for all values of the arguments and k). 

An important feature of Assumption 12.11 is that no control on the 
sizes of r", a^, 6^, and c" is imposed (however, remember that from 
the very beginning they are assumed to be bounded). 

Assumption 2.2. For any unit / G M"' and r/ > 0, we have 

IVa^l <Ko(v^ + r/). 

Remark 2.1. It is easy to see that Assumption 12.21 is satisfied (with, 
perhaps, different Kq) if and only if o"^ := v^a^ satisfies the Lipschitz 
condition with constant Kq with respect to x. 

Indeed, the necessity follows after letting r] I and the sufficiency is 
a direct consequence of the formula 

5,K = 2cr,%,/a," + r/(5,,,cr,")^ 

Below we are also using the well-known fact that a continuous func- 
tion v{x) is Lipschitz continuous with constant K if and only if its 
generalized gradient = D^v satisfies If^^l < K (a.e.). 

Definition 2.2. Let i? be a finite subset of and p{x, y) a real-valued 
function on M.'^ x R'^. For an xq G M.'^ we say that the operator 

S : u ^ Su, Su{x) = p{x, X + y)u{y) (2.4) 

2/esu{0} 

respects the maximum principle at Xq relative to B if, for any function 
0(x) such that 0(a;o + y) > (p{xo) for all y E B, we have Su{xo) > 0. 

Obviously, the operators 6n,i and ^ respect the maximum principle 
at any point relative to appropriate sets. 
For h > set 

Liu = al/\h/,u + hl5h,,,u - c"m, (2.5) 

where and throughout the paper the summation convention is enforced. 
For each t the operator = L^(t, x) can be considered as an operator 
on functions defined on W^. 



6 N.V. KRYLOV 

Assumption 2.3. We have > ho{b'^)~ (recall (l23il ). 

Remark 2.3. It is easy to see that Assumption 12.31 implies that for 
h e (0, ho], t E M, and a E A the operator L'^{t, x) + d^it, x) respects 
the maximum principle at any point Xq relative to Aq, where 

Ao:={/i4:fc = ±l,...,±rfi}. (2.6) 

In turn, provided that all ^k are different, the said property of Lf^{t, x) + 
c°'{t,x) implies what is required in Assumption 12. 3[ 

To satisfy Assumption 12.31 it is sufficient to require that 6^ > 0, in 
which case what we use is just an upwind discretization of the "trans- 
portation" term. 

Remark 2.4. The operators are natural approximations of the op- 
erator 

L^u = ateieiu,.,, + hie^u,. - cu (2.7) 

in the sense that L'^u — > as /i J, for all smooth u. 

One may wonder how wide is the class of operators given in the usual 
form 

Lu = a^^u^i^j + Vu^i (2.8) 

which admit such a special approximation. We discuss this issue in 
Section [9l 

Next, we describe the free term in the equation, which are given by 
a real-valued function 

r = r(p,^,t,x) 

defined on A x M^^i x M x M x R^. 

Assumption 2.4. The function /° is bounded, /" is continuous in a, 
continuous in {p,ip,x) and, for any a and t, its generalized gradients 
Dpf"', D^f"', and D^f"' in p,ip, and x, respectively, satisfy 

\Dp.n<Ko^, k = ±l,...,±di, \D^n<Ko, \DJ''\<K3 

for almost all {p,ip,x) e M^^^^ x R x M'^. 

For fixed /i, r > we consider the equation 

sup[r°5^M + L> + = 0, (2.9) 

where 

9°' = g'^it, x) = P{5h/^u{t, x),u{t, x),t, x). 
Observe that equation (12. 9p takes the form (12. ip . The presence of 
in these equations allows us to treat the normalized Bellman equations 
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(see |9]), which arise, for instance, in optimal stopping problems or 
problems with singular control. 

Fix a vector / G M°' with |Z| < Kq and a number rj G (0, h]. Set 

hi = h for \i\ = l, di, hd+i = /^-(d+i) = ??, £^+1 = = /, 

A = {M±i,...,/id+i^±K+i)}. (2.10) 
We treat Aq as a list rather than the set with specified elements, even 
if £i = £2 we include in the list this vector twice. 

Observe that in (12.51) only ^ are involved. However, the method 
of "shaking" the coefficients requires estimates of difference derivatives 
in all directions and not only along the mesh. This is the reason why 
we introduce A. Set 

n 

A„ = ^ A = {x : X = /i + ... + /„, /i, G A}, Aoo = |JA„ 

1 n 

Mt = {(nr) A T : n = 0, 1, ...} X Aoo, Mr = Mr n {[0,T) x W^). 

(2.11) 

Fix a finite set Q C Air, assume that 

Qlo := g n ({0} X M^) 7^ (2.12) 

and define 

Q = QU{{t + rr(t), x) : (t, x)eQ,t + Trit) = T}, 
Q1 = {{t,x) eQ :t <T,{t + TTit),x) G g, (t, X + A) C g}, (2.13) 

9ig = g\g?. 

Obviously, it may happen that Q = Q. The subscript 1 is used above 
because later on we will need a "fatter" boundary d2Q. 

Finally, define T' as the least nr, n = 1,2,..., such that nr > T, 
recall that m G M (see (12. 3p ) is a given fixed constant and introduce 

m = e"'\ t<T, m=e'^^\ e(+) = evl, e(-)=eAl, 

Cm = , A = inf [c"(t, x) + r"(t)c„] . (2.14) 

7" a,t,x 

Introducing a discontinuous function ^(t) may look unnatural. How- 
ever, what is important for us is that 

^6^u = e-'^-6^{^u)-cU^u) (2.15) 

on Air for any u = u{t,x). 

Everywhere below in this section m is a given function on A^t sat- 
isfying (12. 9p in g. In our first result no control on the sizes of r°, a^, 
6^, and c" is imposed. 
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Theorem 2.5. Let h G (0, /iq]- Then, under Assumptions \2.1\ through 
l^there are constants N = N{di, Kq), N* = N*{di, Kq, K3) such that 
if X > N then on Q\q 

Kiu\ < iV*e'"^(^+-) [1 + max \^^_^u\ + max{\^^.)5,,eM + ■ 

(2.16) 

We prove this theorem in Section |H 

Remark 2.6. This theorem is similar to Theorem 5.2 of [I3] and entails 
all the consequences derived from the latter in [13] and [7j. In partic- 
ular, by using Theorem 5.6 of [13] and comparing the equations for u 
and u{t+ -, ■) an estimate of 5rU can be obtained if we require the data 
to have bounded derivatives in t. 

Remark 2.7. Theorem 12.51 has an immediate application to elliptic 
equations. In that case u is independent of t, one can take r° = 0, 
and use as large negative m as one wishes without affecting A. Then 
it is seen that in (12.161) the maximums over Q and diQ reduce to the 
maximums over Q\o and Q\q fl diQ, respectively. 

Our second result is about Bellman equations with more general 
quasilinear operators. This time (12. 9p is assumed to be uniformly non- 
degenerate in the space generated by f^'s. We will allow f°'{p,ip,t,x) 
to grow quadratically with respect to p and therefore no 6^ are needed. 
The term also could be absorbed in However, we keep it, in 
order to state Theorem 12.111 in a simpler way. 

Assumption 2.5. (i) The functions also depend on if): 

and equation (12. ip holds in Q, where F is defined by (12.20 with a'^{ilj,t,x) 
in place of a^(t,x). The functions a'^{ip,t,x) are Lipschitz continu- 
ous in X with constant K3, Lipschitz continuous in ip with a constant 
UJ G (0, 00), 

flfc > ^, 1^1 < di, > -K3. 

(ii) The function is continuous in a, continuous in {p,ip,x), and 
for all values of the arguments, satisfying < Ki and \p\ > K2, it 
holds that 

\f^\<u;\p\' + K,. 

(iii) For each a and t the generalized gradients Dpf°', D^f", and 
D^f"' of /" with respect to p, ip, and x, respectively, satisfy 

W\ < uj\p\ + Ks, |/^^r I < u\p\^ + Ks, 

\DJ"\<u\pf + K, (2.17) 



A PRIORI ESTIMATES 



9 



(a.e.) on the set {\ip\ < Ki}. 

Remark 2.8. Clearly, Assumption 12.51 (ii) is satisfied with any uj > 
and appropriate K^^u) if 

sup \r{p,^,t,x)\ = o{\p\^) 

as \p\ oo. This includes all functions affine in p provided that 
the coefficients are bounded. Similar situation occurs with Assump- 
tion O (iii). 

Assumption 2.6. For a constant C > 4 depending only on di, the 
exact value of which can be determined by examining the proof of 
Theorem 12.91 we have 

CKi{l + Ki)u <6. (2.18) 

Theorem 2.9. Let 6^ = and let Assumptions \2. 1\ \2.^ and \2.6\ be 

satisfied. Assume that \u\ < Ki in Q and \6h/t.u\ < on diQ if 
\k\ < di. Then in Q 

<N = N{d^, 6, K,, K^, Ks), \k\ < dr. 

In particular, N is independent of T . 

The proof of this theorem is given in Section [51 

Remark 2.10. If uj is large, we need Ki to be small in order to satisfy 
fl2.18p . that is, we need u to be small. By replacing u with u — •y, 
where 7 is any constant, we see that, actually, we need the oscillation 
of u rather than u itself to be small if u is not. This restriction could 
be completely avoided if we proved an interior version of Theorem 12.91 
and a priori Holder continuity of u. It seems to the author that this is 
possible, but requires much more work. 

Theorem 2.11. Under the assumptions of Theorem \2. ^ suppose that 
\^r),iu\ < on diQ and are independent of ip. Then there is a 
constant N = N{di, 6, Ki,K2, K^), such that if X > N, then on Q\o 

This is a simple corollary of Theorems 12.91 and 12.51 with Hq = h in 
the latter. Indeed, once we know that the values of \6h/kU\ and |m| 
are dominated by a constant, the behavior of f°'{p.,ip.,t^x) for large |p| 
becomes irrelevant and we can even multiply it by an appropriate cut- 
off function in such a way that the new would satisfy Assumption 
12.41 and u would still satisfy the new equation. 

Our next result is about second- difference estimates. 
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Assumption 2.7. (i) The function is independent of p and ip. 

(ii) For any i,j = ±1, ...,±{di + 1) and standing for any of the 
functions 6^, c°, and we have 

(2.19) 

Typical case when the third inequahty in (12.191) is satisfied occurs if 
a1 = (o"^)^, where cr^ is bounded and twice continuously differentiable. 

In contrast with the above results in which no control on the magni- 
tudes of r", a^, 6^, and is required, this time we need the following. 

Assumption 2.8. We have 

6<snpat<Ko, r^,\bt\Ac%\n < K^. 

The following assumption is about a special structure of the set of 
our basic vectors ik,k = ±1, ...,±di. For di = 2 and the standard 
grid (generated by ±ei,±e2) it means that this set contains all eight 
neighboring points of the origin on the grid. 

Assumption 2.9. There exists an integer 1 < do < di such that for 
the list 

C:={hi±i,...,hi±do} (2.20) 
and any ik with di < \k\ < di there exist li,l2 & C such that 

One may think that Assumption 12.91 excludes the equations with 
only one spatial variable, where it is natural to take di = 1 and Aq = 
— However, we do not require Ik to be nonzero, and one can 
take Ao to be {^1,-^1,4,-4,^3,-4} with 4 = and 4 = 4- In 
that case Assumption 12.91 is satisfied with £ = {£i, — £1, 4, — 4}- By 
the way the fact that now the origin is one of Ik in no way contradicts 
Assumption 12. 8[ because in that case 5h,i^(t) = and one can assign any 
value to without changing the equation. 

Define 

Ql = {{t,x) eQ ■.t<T,{t + Trit), x) e Q, {t, x + Aq + Aq) C Q}, 

d2Q = Q\Q°2- 
Here, naturally, x + Aq + Aq = {x + y + z : y, z & Aq}. 

Theorem 2.12. Suppose that Assumptions l2Jl {MM IMIM W^W^ 

satisfied. Then there exists a constant N = N{6,di, Kq) such that if 
X > N, then in Q\q for i,j = ±1, ±di we have 

\5h,iA,eM < N^e'^^^'^-^-^R, (2.21) 
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where N* = N*{ho, 6, di, Kq, Ks), 

R = 1 + max + max(^(_)(5j'u)~ 

+ max max _)6h,e,Sh/M + max max \^/.)6h,iM)- 
\i\,\j\<di d2Q ' ^ ' " ' \i\<di Q ' ^ ' 

This theorem is proved in Section [7] following a quite long Section [6] 
that contains the proof of Theorem 12.121 under additional assumptions. 

Remark 2.13. To get "closed" estimates of Sh/^Sh/iU we need to exclude 
Sh/iU and S^u from R. This can be done by using Theorem 12.51 and 
the idea from Remark 12.61 Another situation when S^u drops out 
presents when u is independent of t, so that, actually, we are dealing 
with elliptic equations. We say more about this in the comments after 
Theorem 12.141 

In case of independent of x Assumption 12.91 is not needed. 

Theorem 2.14. Suppose that Assumptions {MJj \23, IMIM andlKE 

are satisfied. Also assume that are independent of x, {Shi/ib'^l < 
^^V^' N ^ 'j^i + 1; 1^1 ^ di, and A > 0. Then in Q\o for k = 
±1, ±di we have 

|A,,4m| < Ar*e'"^(^+^)i?o, {\,iu)- < N*e^^^^^^'^R, (2.22) 

where N* = N*{\, ho, S, di, K3), 

Rn = 1 + max \C(-)u\ + max(f(^)6^u)~ 

+ max max \^/)Ah^^iM\ + max max \^,_)6h„eM)^ 

\i\<di diQ \i\<di Q 

and R is obtained from Rq by taking di + 1 in place of di. 

This theorem proved in Section [S] is a direct generalization of the 
corresponding result from |5j: lower order coefficients are allowed to 
depend on (t, x) and we consider parabolic equations. In connection 
with the latter observe that if r"" = (elliptic case), then one can let 
m — s> — 00 and see that in the definitions of Rq and R one can replace 
diQ with g|o n diQ. 

3. Some technical tools 
For any rj G M*^, u > set 

T^,^ip{x) ■=i){x + vif]). (3.1) 
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Lemma 3.1. For any u > 0, li,l2 G K'^, and functions a{x),ip{x) 

+ h[5y^i^a + 5y^i^a]5u,i25u,i^^l) + {5y^i^5y^i^a)Th,i^+i2^, (3.2) 
In particular, 

This lemma is proved by straightforward computations (cf. [T3]). In 
the following lemma we use Definition I2.2[ 

Lemma 3.2. If an operator 

Sil){x) = p{x,x + y)il){x + y) 

j/eAu{o} 

respects the maximum principle at a point xq € relative to A anc? 
if) is a function such that ip{xo) < 0, then —Sip < S{iIj~) at xq- In 
particular, (j)~S(j)~ > —(j)~S(j) at xq for any function (p. 

This follows from the definition and the fact that ip + tp~ > on 
xq + A and + ^Z'" = at xq. 

The following lemma from [13] is used in the proof of Theorem 12.121 

Lemma 3.3. Let ip be a function on M^, u > 0. Then 

|A,,^^| < + (3-3) 

4. Proof of Theorem 12.51 

We start with some preparations. From now on index k will run 
through {±1, ±di} and i,j through {±1, ±{di + 1)}. By and 
N* in this section we denote generic constants depending on the data 
as in the statement of the theorem. We use the notation fl2.1ip through 
(I2.14P and introduce few new objects. We need two constants e and n 
defined by 

e~^-2edi = l, 4/i = (di + 1)"^ A e. 

Introduce 

r = {7 = (7. : ^ = ±1, ±K + 1)) : 7. e [e, e"!]}, 
By using fl2.15l) we see that in Q 
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Also introduce 

= P^v, Vi = 6iV, 

i 

Observe that 

Finally, let (70,^053:0) G F x Q be a point at which V^^ attains its 
maximum value over T x Q. 

Theorem 4.1. The assertions of Theorem \2.5\ hold true if in addition 
to its assumptions (to,a;o) G Ql and 

v-{to,xo)>{l/2)mjix\vi\. (4.1) 

To prove this theorem we need an auxiliary result. 

Lemma 4.2. Assume (14.11) . Then the operator P^^^^ respects the maxi- 
mum principle at {to, Xq), that is, for any function (p such that </)(xo) > 
(j){xo + rj) for all rj E A, we have PyQ^(f){tQ, Xq) < 0. 

Proof. Since P-y/^l = we may assume that 0(xo) = 0. Then at 

Piof^(p = (%7oi - fJ^Vi)(j){xo + hdi), 
which is negative since (f){xo+hiii) < and v~^'yoi—fiVi > f ~ (e— 2/i) > 
(A'- ^ ^/2)- The lemma is proved. 

We also need the following construction. Notice that, if (to? 2:0) G Q, 
there is a sequence an E A such that at (to,xo) 

lim [r""(5^M + a^" A^m + b'^"6kU - c""u + g°"] 

n— >oo 

= sup[r"5^M + a°AkU + b'^dkU - c'^u + = 0. 

Since the numbers of possible values of t for points in Q is finite, and 
the functions (t, x), b^{t, x), c"(t, x), /"(p, V'? t) ^) are uniformly con- 
tinuous functions of {p,ip,x), there is a subsequence {n'} C {1,2,...} 
and functions f, ak(t,x), bkit,x), c{t,x), f{p,ip,t,x) such that they 
satisfy our assumptions changed in an obvious way and 

(f, afe(t, x),bkit, x), c{t, x), f{p, ip, t, x)) 

= hm (r°"'(t),a°"'(t,x),6^"'(t,a:),c°"'(t,x),r"'(p,^,t,x)) 

n'— >oo 

on Q for all p, ip. 
Obviously, for 

g{t,x) := f{5ku{t,x),u{t,x),t,x) 
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at (to, xq) we have 

fS^u + dk^kU + bkSkU — cu + g = 0, 

e-'^^fd^v + akAkV + hdkV - (c + fc„^)v + = (4.2) 

and, if (tQ,Xo) G Q1, then for any i (= ±1, ...,±{di + 1), the shift 
operator T is introduced in (13.11) ) 

Th„e, [e-^^^fd^v + dk^kV + hdkV - (c + fCm)v + < 0, (4.3) 

where and below for simphcity of notation we drop (to,a^o) in the ar- 
guments of functions we are deahng with. 
Proof of Theorem 14. 1[ Set 

Ll = flfcAfc + hSk, Lh = Ll-{c + fCm)- 

By Assumption 12.31 Lemma [XT] and Lemma IX^ (with in place of 0) 

> = 2v~Llv- + 2^iViLlvi + h+ fih 



where 



Ji := 2ak{SkV^^y + /^fefe(4%)^ 
I2 := 2dk ^{SkVif + hbk ^{hvif. 

i i 

Since 2dk + hbk > dk we conclude 

%^Xo - f^^iL^Vi > il/2)fidk J2(^kViY- (4-4) 

i 

On the other hand, by (14. 2 p and (14.31) and Lemma [4.21 at {to,xo) 

P^,^[e-"'^f6^v + LhV + ig] < 0. (4.5) 
Owing to (14.41) we obtain 

P^of^LhV = t;:^jI°t;^o - mllvi + (c + fCm)Vy„^ + 13 + ^ + 15 

> XV^,f, + {l/2)l2dk J2(^kVi)^ + h + h + h: 

i 

where 

/g := v~^{P.,jLk)i:^k'v + hiV~^^-^Qi{bidk)IAkVi 
-l2Vi{6idk)AkV - j2hiVi{6idk)AkVi, 
h ■= [%^o» ~ fJ'Vi]iSibk)Th„£,Vk, 

Upon observing that 

Afc = -6k6_k, hAk = 6k + 6_k, hi < h 
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and by assumption (14. ip 

h\6kV_k\ < 2m8ix\vi\ < 4f_, \h^AkVi\ < 4max|fj| < 8v', 

Q,i Q 

we find 

l^sl < Nv::^^{y/U^ + h)\SkV_k\ + y^.iV^ + h)\hAkVi\ 

i 

< Nv- J2 V^k\SkV^\ + N{v-J' < N{v~f + (l/4)/xa, J2(6kV,f. 

i i 

This yields 

P^ot^LhV > XV^,^ - N{v-J^ + (l/4)/iafc J^^SkV,)^ + h + h. 

i 

Next, obviously, 

I/4I < N{v~^)\ \h\ < N*v- max \v\. 

Q 

Therefore, and since (v^^)^ < ^om' 

Pjof^L^v > (A - N)V^,^ - N* max \v\^ + {l/4)fiak J^i^kVif- (4-6) 

Now we deal with other terms in (14. 5p . Note that, since (to + 
rT{to),xo) G Q, 



so that (cf. Lemma 13.21) 

By recalling that r° > 0, we find from (14. 6 p and (14.50 that 

N* max \v\^ > (A - N)V^,^ + iP^.^g + (l/4)/^afc ^^(^fet;.)^ (4.7) 

By Assumption 12.41 

i 

It follows that the sum of the last two terms in (14. 7p is greater than 
-NV,,, - N*iV^/^^ > -NV,,, - N*e 

Hence 



iV*(max + > (A - Ni)V 
Q 



70 A*' 
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which for X — Ni > 1 shows that on Q\q 

< iV*(max \v\^ + e) < iV*e('+)(T)(max + 1). 

This imphes fl2.16p and the theorem is proved. 

In hght of this theorem to prove Theorem 12.51 we only need to show 
that in case the assumption of Theorem 14.11 are not satisfied one can 
obtain the assertion of Theorem 12.51 differently. To do that we need 
two lemmas. 

Lemma 4.3. Take a function on M'^ and assume that 

max < max|(5i0(O)| 

i i 

for X G A. Then 

max|5i0(O)| < max ({P^(p{x)y : 7 G T, a; G {0} U A) . (4.8) 

i 

Proof. Take a j such that 

max|5i0(O)| = \5j(j){0)\ 

i 

and first assume that 

15,0(0)1 = -5,0(0). (4.9) 

Then take 7,- = and 7^ = £ for i 7^ j. Since — 5,0(0) > 5j0(O) and 
— 5,0(0) > 0, we have 

P^0(O) = £-^5,0(0) + eJ2 

< {e-^ - 2edi)5,-0(O) = -(^"^ - 2€di) max |5i0(O)| 

i 

and (14.81) follows since — 2edi = 1. 
If ( 14.91) is not satisfied, then 

max|5i0(/i,£,-)| < |5,-0(O)| = 5,0(0) = -6^j(f){hjij), 

i 

which combined with an obvious inequality between the extreme terms 
yields 

max 1 5j0 (/;,,•£,•) I = —6^j(t){hjij) = |5„,-0(/i,-£j)| 

i 

= |5,-0(O)| = max|5i0(O)|. 

i 

By the above argument applied to the point hjij in place of 
minP^0(/i,£,) < (e^^ - 2edi)6.j(f){hjij) 

= -{e-^ -2£rfi)max|5i0(O)| = - max |5i0(O)|. 
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The lemma is proved. 

Lemma 4.4. Condition (14.11) is satisfied if {to,xo) G Q" and 

maxvf = ma.xv^. (4-10) 

Q,i Qi^i 

Indeed if (14.1 p does not hold, then 

max[t;~]^ < V^f^f,(to,xo) 



< (2(rfi + l)fi+ 1/4) max = (2(rfi + l)/i + 1/4) maxt;? 

which contradicts Lemma [4.31 since fx < l/{4{di + 1)). 
Proof of Theorem 12.51 If {to,xo) E diQ, then in Q\o 

\6,,u\ = \6,,v\ < fx~'/X£ < f^-'^Xl"ito,xo), 

where the last term is obviously less than the right-hand side of (I2.16p . 
Furthermore, if (to,Xo) G Q" but (14. ip is violated, then by Lemma l44l 

maxu^ = maxf,^ 

Q,i diQ,i 

and one can use the above argument. Finally, in the remaining case 
Theorem 14.11 is applicable. The theorem is proved. 

Remark 4.5. As in ^ one can relax conditions on A by adding in its 
expression any large constant times min^ a^. This can be shown by 
considering maximum points of + z/f ^, where z/ is a large constant. 
This would also allow us to relax the condition on D^f to \D^f\ < 



5. Proof of Theorem 12.91 

Our goal is to show how to choose an appropriate C = C{di) in 
(I2.18p . Below in this section by N we denote generic constants de- 
pending only on di, 6, Ki,K2, but not oj. 

First of all, observe that / does not enter either equation (12. 9p or the 
statement of the theorem. It is involved, however, in the definition of 
diQ making it "fatter" . Because of that if we additionally assume that 
/ = 0, the result will be stronger. Therefore, we assume that / = 0. 

Set m = and introduce e, F, f^, fj, VC^^, P-y, P-y^ as in 
the beginning of Section |H However, since v = u, we also write 
Uy, Ui, and U^^ instead v^, Vi, and V^^, respectively. Since / = 0, 
■^±(^1+1) = "^itCdi+i) = and now there is no need to allow i to take the 
values ±{di + 1). Therefore, we restrict it to the range ±1, ±di. As 
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everywhere in the article index k runs through ±1, ±di. This time 

we take 

8/i = d^^ A e. 
Introduce k > as a solution of 

max[f/^^ + Ku^] = AkKI. (5.1) 
rxQ 

Observe that if k = is a solution of (15. ip . then = and the 
assertion of the theorem is trivial. Therefore, without losing generality 
we may assume that for k = the left-hand side of (15.11) is strictly 
greater than its right-hand side. Furthermore, as a function of k the 
left-hand side is convex increasing with Lipschitz constant not greater 
than K\. It follows that (15.11) has a unique solution k > 0. After that 
we define (70, t^, Xq) as a point in F x Q at which the maximum in (15.11) 
is attained. For simplicity of notation we drop the arguments (to, Xq) in 
what follows. We also use the abbreviated notation 5i, introduced 
in Section m 

First we show that [m^J^ is the main term in U^^^^ + nv?. 
Lemma 5.1. Assume that (14.101) holds. Then 

[u-^f > 2kKI (5.2) 
2die~^ Yii&5c\ui\ > > (1/2) niax (5-3) 

i ^° Q,i 

In particular < and the operator P^^^ respects the maximum 
principle at (tg, Xq). 

Proof. First, notice that, owing to (14.101) . Lemma [4.31 and the defi- 
nitions of K and {tQ^Xo) we have 

AkKI = + /iVw^ + km^ < [mZj]^ + 2c/i/im_axM^2 + kK^ 

Q,i 

I 

< [u-f + 2rfi/imax[Mr]2 + kK^ < [u-f + Id^^iAnKl + KKj. 

' rxQ ' ' 

This implies (15.21) since Sdifi < 1. 

The first inequality in (15. 3p is obvious. If the second one is wrong, 
then 

max[n^]^ < U-y^^n + f^u^ < [2difj. + 1/4] rnaxw^ + kK^ 

TxQ Q,i 

< [2difi + 1/4] max[M-]2 + kK^ < (1/2) maxiwl]^ + nK^ 

VxQ ' FxQ ' 

contrary to (15. 2p . The last assertion of the lemma follows from Lemma 
?2l The lemma is proved. 
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Next, if {to,xo) G diQ, then 

AKKf = Uy^f, + Ku^ <N' + kKI 3KKf < N', 

and u\{t,x) < (4/3)yU~^A^' in Q according to (15. ip . In this case the 
assertion of the theorem is true. Similarly, as in the proof of Theo- 
rem 12.51 we get the result if (I4.10p is violated. This justifies the first 
two assumptions in the following set which we impose: 

ito^xo) e Ql, rnaxM? = maxw?, 

Q,k Q1,k 

Adie'^W^/^ > max ImJ, W > + 8X1X36'^ + 1. (5.4) 

The third relation in (15. 4p follows from Lemma [5?T] and the last assump- 
tion in (15. 4p restricts us to the only nontrivial case in light of (15. 3p . 

Next again owing to the fact that the number of points in Q is 
finite we can find some functions f{t), (ikit, x), c{t, x), and f{p, ip, r, t, x) 
satisfying Assumptions 12.11 and 12.51 and such that for 

g{t, x) := f{uk{t, x), u{t, x), t, x) - c{t, x)u{t, x), Lh := a^Afc 

we have (recall that the arguments (to,xo) are dropped) 

f6^u + LhU + g = 0, Th/^ [f6^u + LhU + g] < (5.5) 

for any i (= ±1, ±di). 
Now, as before 

> Lh[U^of^ + Ku^] > -2u~^LhU^g + 2fiUiLhUi + 2KuLhU 

+ 2afc[(4n~ )^ + /i ^{6kUif + K{6kuf]. (5.6) 

i 

Furthermore, 

i 

> -2u~^5'^u^^ + 2fiUi6'^Ui + 2ku5'^u. (5.7) 

We multiply (15.70 by f, add the result to (15. 6p . and use that f > 0, 
c > — -R's, flfc > 5, |m| < Ki, and 

u{f5^u + Lhu) = —ug = cu^ — u{g + cu) > —K^Ki — u{g + cm), 

where |(7 + cu| < ujW + (recall (I5.4p ). Then we obtain 

u{f5^u + Lhu) > -ujKiW - 2K1K3, 

u- [f(5^ + Lh]u^^ - i^Ui[f6^ + Lh]ui > 6fiM + kW{6 - K^uj) - 2KK1K3, 
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where 

M := Y,^5uUi)\ 

i,k 

Since C > 4 in Assumption 12. 6[ we have Kiu < 5/4 which along 
with ( 15I1) leads to 

W{6 - Kiio) - 2K1K3 > {3/4:)SW - 2K1K3 > {6/2)W, 

[f(5^ + Lh]u^, - fxUi[fS^ + Lh]u^ > 5fxM + {S/2)kW. (5.8) 
On the other hand, owing to (15.41) . (15. 5p . and Lemma [5. II 

P^,^[f5ju + LhU + g]<0. (5.9) 

Here due to ([EH]) 

P^f,f,[f6ju + Lhu] = u:^^ [f6^ + - fiUi[f6j + Lh]ui + I 

> 6ijM +{1/2)6kW + I, 

where 

/ := [u:^g7oi - fiUi]{6iak)[AkU + hAkUi]. 

Below by C we denote generic constants depending only on di. It 
follows from the estimates 

\6idk\ < + uj\ui\, \hAkUi\ = \6kUi + 6-kUi\ < 2M^^'^, 
\Aku\ < M^'\ |m;„7o^ - < CW^'\ 

that 

|/| < CW^/^M^/'^^K^ + uW^'"^) < NW + (l/2)5/iM + Cd'^uj^W^. 
Hence, 

P^,^[fS^u+Lhu] > {l/2)SfiM+{l/2)6KW-NW-C6-^ujV. (5.10) 
To estimate P'y^fj.g recall that c > — A^, U^^^ > and observe that 
-P-yotiicu) = -cPj,,f,u - [u:^^^-fok - ^Uk]{Th^i^u)5kC 

= cU^,^ - [m- 7ofc - fiUk]{n,e,u)6kC > -NW - NW'^^ > -NW, 
where the last inequality follows from (15. 4p . Furthermore, 

h6i{g + cu) = f{Th/Mk, Th^e.u, to, xo + Mi) - f{uk, u, to, xo). 

Owing to Assumption 12. 5[ (15. 4p . and the mean value theorem (this is 
the place, where one cannot assume that (12.171) holds only for large p) 

Mg + cu)\< CM'/\ujW'/^ + Ks) 

+C[W^^^{ujW + Ks) + ujW^^^ + K3]. 
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Note that the coefficients of 6i{g + cu) in Pj^^f^^g + cu) are dominated 
by W^^'^ and for any p > 

uM^/^W <pM + p-^uj^W\ KsM^/^W^/^ < pM + p-^NW. 

Therefore, 

\P^,U9 + cu)\< {1/2)6 fiM + C[rV + uj]W^ + NW, 

> -{l/2)6pM - C[6-^uj^ + uj]W^ - NW 

and dSl]) and (ISTTOj) yield 

{1/2)6kW < C[6-^uj'^ + uj]W^ + NW, 

26k<Ci[6-^uj^ + uj]W + Ni. (5.11) 

Since W < p^^AKfn (see (15. ip ). we see that if C in Assumption 12.61 is 
such that 

C > 4p-^Ci, 
then (recall that C > 4 and uKi < 6) 

Ci[5-^uj^ + uj]W < CKf[6-^iu^ + iu]K< C[KiUJ + KIuj]k < 5k. 

In this case (15. lip allows us to conclude that k < 6~^Ni and we get 
the assertion of the theorem from (15. ip . The theorem is proved. 

6. Conditional estimates of the second-order differences 

In this section we suppose that all the assumptions of Theorem 12.121 
are satisfied apart from Assumptions 12.81 and 12.91 The notation in this 
section are somewhat different from Sections H] and [51 Of course, we 
use our basic notation from Section [21 for instance, ^ and A are defined 
in dsn. 

For e G (0, 1] set (observe that now ■jk = 7__fc) 
T{e) = {7 = (7fc : = ±l,...,±di) : 7^ = 7.^,5 < 7/, < e~\Vfc}. 
Fix a constant p = p{di,e) > such that 

Wdjp < e\ 

In this section the indices i,j,k,p,q run through {±1, ...,±di}. The 
main result of this section is the following. 

Theorem 6.1. Assume that 

3/xmax^ [(C^m/v.w)^^ < max [{^^-f,Ah,e^u)'Y . (6.1) 

Then there exists a constant N = N{6, p, di, Kq) such that if X > N, 
then in Q\o for i,j = ±1, ±di we have 

\Sh,e,SH,iM<N*e'^^^^+^\l+ max \^^^)6h,e,6H,i,u\ 

d2Q,p,q 
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+ max + max \C,(_)6h,epU\) , (6.2) 

Q Q,p 

where N* = N*{hQ, e, /x, di, Kq, K^). 

Below in the section by and A^* we denote generic constants of 
the same type as in the theorem. As before, we use the abbreviated 
notation 

Introduce v = as in Section H] and fix a constant z/ > 1. Set 
= 7iAj, = P^v, Vi = 6iV, Vij = 6j6iV, 

P^t,u4> = v:^P.,(f) + fiv~j6j6i(j) - uviSi(f), 

Observe that this time again P^^uV = —V^^u and also note that (16.11) 
is equivalent to the following 

3/imaxW2 < max {v~Y. (6.3) 

We introduce (7o,io;2;o) as a point in V{e) x Q maximizing V^^y 
and first prove few auxiliary results. Below, as usual, we drop the 
arguments {tQ.Xo). 

Lemma 6.2. (i) For {t,x) G Q" and any i,j 

\vij{t,x)\ <maxW2^'^. (6.4) 

(ii) If (16. Sp holds and {tQ,Xo) G Q and 

i/ max M/^i < /X max 14^2, (6.5) 

then at (to, Xq) 

fimaxW2 < [%]^; umaxWi < (6.6) 
Furthermore, if additionally, 

h\fv < £, 

then the operator P-^^^u respects the maximum principle at [to, Xq) rela- 
tive to Aq + Aq, that is, for any function such that 4>{xo) > 4>{xo + rj) 
for all rj E Aq + Aq, we have P^^^u4>{'to, xq) < 0. 
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1 II 

Proof, (i) Obviously fZj^- < W2 on Q. Since (no summation in 1) 

Th,e,vZij = vfj, we get vfj{t,x) < W\^'^{t,x^ Mi). This proves (i). 

(ii) The second estimate in (16. 6p follows from the first one and (16. 5p . 
Assuming that the first estimate in (16. 6p does not hold, we obtain at 

(^o,a;o) 

'V-ynu.v < 2u max 14^2 + z/maxVTi < 3u max 14^2, 

' Q Q Q 

contrary to (16.30 . This proves (16. 6p . 

To prove the last assertion of the lemma we take a function with 
described properties and without loss of generality assume that 0(xo) = 
0. We also note that 

h'^6i6j(j){xo) = 0(a;o + hii + hij) - 0(xo + hii) - (j){xo + hij) 

+4>ixo) < -0(a;o + hii) - 0(2^0 + h£j) 
and Aj0(xo) < 0. Therefore, as usual dropping the arguments {tQ,xo) 
in V,,,, we infer from (16. 6p that 

h^P-yoiiAih, Xq) < iq^E ^(0(xo + hii) + 0(^0 - hii)) 

i 

-HV~j{(f){xo + hii) + 4>{xq + hij)) - hvvi(j){xQ + hii) 

i 

-v^l](0(a;o + /i^i) + 0(a;o + /i^i)) - /iv^^0(xo + hi,)] 

i 

0(xo + hii)[2e -2di^- h^. 



^0 



The last expression is less than zero in light of the fact that 4(ii^//i < e 
and h^Jl' < e. The lemma is proved. 

Remark 6.3. This lemma can be generalized to the case when ^^'s come 
with different /ife's, but the h^s should be comparable. This is the 
reason why in Theorem 12.121 we do not include h^^^i and tdi+i- 

Set 

Zk := v:^^AkV^^ + fxv:r.AkVij - uViA^Vi, 
Zk := %4^^7o + f^v~j6kVij - uViVki, 

id i 

These objects evaluated at (to, Xq) will be extensively used below in the 
section. 



24 N.V. KRYLOV 

Lemma 6.4. // {tQ,xo) G Q", then at {to,xo) we have for any k 

>Rl + Rl, + /ii?^ + ^^Rl, + vRl + vR% > 0. (6.7) 
Furthermore, for any a & A and h < /io/2 

4«Z, + b^Zk) > 2alZ, + al{Rl + /xi?^ + vRl). (6.8) 

Proof. The first inequality in 06.71) follows from Lemma 13.21 To 
prove the second one it suffices to observe that 

> ^kV^af^u = 2v~^^kV~^ + 2^Vlj^kVij + 2l'Vi/^kVi 

+Rl + R\ + /ii?^ + ^iR^_^ + vRl + vR%. 
Next, using that h < /io/2, by (16. 7p and Assumption 12.31 we get 
> «Afc + 2hl5k)V,,,, = 2v-[atAk + 26^5,]i;- 

+2fiv7^[atAk + 2hl5k]vi^ + 2pv,[al/\k + 26^4]t;, 
+2a^[/?^ + /ii?^ + uRl] + 2/i6^[i?^ + /ii?^ + uRl] 

>-2v-^[al/\k + 2hl5k]v,, 
-2iivr.[al/\k + '^hl5k]vij + 2i/t;iKAfc + 26^4]t;, 

+2at[Rl + /"^fc + ^^^.1 + '^hhl[Rl + /^^^ + ^RV\ 
> -2alZk - Ahlzk + al[Rl + /ii?^ + vR^] 

and (16.81) follows. The lemma is proved. 

In the following lemma we do the most important step in the proof 
of Theorem 16.11 Set 

Wi := maxWi. 

Q 

Lemma 6.5. Under the assumptions of Theorem \ 6.1\ there are con- 
stants N, N*, 

V = i^*{ho, e, /i, di, Kq, Ks) > 1, h* = h*{ho, e, /i, di, Kq, K^) > 

such that h* < ho/ 2 and, if condition (16. 5p is satisfied and /i G (0, h*] 
and (to,xo) G Q2, then at (to,a^o) for any a & A we have 

J := P,,,,{atAkV + hl5kv) > -N[irf - N*W,. (6.9) 
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Proof. We fix an a G A and drop the superscript a for convenience. 
By f l3.2p (no summation in k,i,j) 

SjSi{bkVk) = h5kVij + {5jhk)vki + {5ibk)vkj 
+h[{Si + 5j)bk\5kVij + {5j5ibk)Th,i^+i^Vk. 
Also by using (13. 2p and the formulas = a_fc and hA^ = Sk + and 
summing with respect to k (but not in i,j) we get 

5j5i{akAkv) = akAkVij + {6jak)AkVi + {6iak)AkVj 

+2[{Sj + Si)ak]6kVij + {6j6iak)Th^e,+eAkV- 
While applying this formula to Aj it is also useful to observe that 

[(5_i + 5i)ak]5kVi_i = -h{Aiak)Ai5kV = -{Aiak){vki + Vk-i). 
Hence, (recall that 7oj = 70 -« and £j = —(--i) 

P^^{akAkV + bk^kv) = -joi6-i6i{akAkV + bkSkv) 
= (afcAfc + bk6k)vy^ - 2-foi{6-iak)AkVi + A-foi{Aiak)vki 
+7oi(Aiafc)Afcf + 2-foi{6ibk)vki + -foi(,Aibk)vk. 
Also everywhere 

v~j5j5i{akAk + bk5k)v = vl^{akAk + bk5k)vij + 2vl^{5jak)AkVi 

[(^i + k)a'k]AkVij + v:r.(^Sj6iak)Th/^+iAkV 
+2vlj{5jbk)vki + 2hvl^{5ibk)6kVij + vlj{5j5ibk)Th,i^+i^Vk. 
Therefore at (to, 2:0) we have 

J = ttkZk + bkZk + I1 + ... + h, 

where 

Ji = -2'joiV~^{S.iak)AkVi + 2^vl^{5jak)AkVi, 
h = fihv'j [{6j + 6i)ak]AkVij, 
h = %[47oi(Ajafc)t;fci + 7oi(Aiafc)Afct;] 

+fiv^j{6j6iak)Th/^+iAkV, 
h = 2v~^joi{6ibk)vki + v~^-foi{Aibk)vk 
+2jjv:r.{5jbk)vki + 2ij,hv~j{5ibk)5kVij + Hv:r.{SjSibk)Th^e,+e,Vk 
-iJVi{5iak){AkV + 2vki) - i^Vi{Sibk)Th/^Vk. 
For h < ho/2 it follows by Lemma [6.41 that 

AJ > 2akZk + ak{Rl + I^Rk + ^^1) + + - + ^h- (6.10) 
Estimating Ii. Note that owing to (16. 6p 

|4/i| < Nv-^{y^+h) \^kVi\ 
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and by (EJ]) 

i i 

< Niv-f + (l/3)/ia,i?^ 

Furthermore, by the formula hAk = Sk + 5h/_k and Lemma 16.21 we 
obtain 

v-^h J2 \^kv^\ < 2v-^ Yl l^'^^l ^ < Niv-^y. 

Thus, 

\Ah\<N{v-y + {l/3)f,akRt. (6.11) 
Estimating 12- Observe that 

I4J2I < NhvTj{,/E;+h)\AkVij\ < hi + /22, 

where (see Lemma [62] and recall that (to^a^o) G Q2) 

k k 

and by the formula \6\ = 6 + 26~, 

J22 = Nhv:r.,y^\AkVij\ = Nhv'^^/a^AkVij 

+2Nhv:[j^/a^{AkVij)~ = Nhy/a^Zk + Nhvl-^/a^{AkVijy + hz, 
where 

J23 = -Nh^/a^{v~^AkV^g - vviAkVi). 
Furthermore, by Lemma [3.21 

Nhv~jy/a^{AkVij)~ < Nhy/a^v~j\Akv:r.\ = Ny/a^v~j\{Sk + (5-fc)^i^l> 

which is majorated by the right-hand side of (16.111) . 
To estimate hs we use Lemma [3.21 to get 

<Niv-f + {l/3)akRl (6.12) 
Furthermore, by assumption ( 16. 5p 

Nhvvi^/a^AkVi = Nvvi^fa^Vki < N{v~^f + {l/?>)uakRl■ 
\t follows that 

/23 < iV(^o)' + {l/3)ak[5kV-^f + (l/3)z/afcK. 

Hence 

I4/2I < N{v~^f + Nh^kZk + (l/3)afc(i?^ + /ii?^ + uRl). (6.13) 
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Estimating I3. We use the following result of simple computations 

Th,e^+i^AkV = -Vk-k + Vkj + v-k,j + Vki + V-k,i + h?AkVij. 

This shows new terms entering J3. All of them apart from the last one 
are similar to the ones which are written explicitly in the definition of 
/a and we show how to estimate only one of them. By Assumption 12.71 
and Lemma [6.21 we have 

\^iv~j{5j5iak)vk-k\ < N{v~^^f+N*\v~^^\ ^i/o^l^fcil < N{v~y+N*akRk- 

k,i 

To estimate the remaining term in I3 we proceed as in estimating I2. 
We have 

^h'^v~j{6j6iak)AkVij 
< Nh\{N + N*^k)\AkV^,\ = hi + /32. 
Here by Lemma W2\ and because (to,Xo) G Q2 

/31 = Nh%\AkV,,\ = Nv~^\{n,e, - 2 + T,,,_Jt;,,| < N{v-J'. 
Next, 

J32 = N*h^v:[jy/a^\AkVij\ = J321 + I322, 

with 

/321 = N*h'^iJ,v~jy/a^AkVij, 
h22 = N*h'v-^^k{AkVij)~ < N*h\^k\AkV-^\ 
= N*hvij^k\SkV-j + S_kV-,\ < N*Wi + {l/3)fiakRt, 
where the last inequality is true since hv~j < 2W^ . Also observe that 

-^'321 = N^h'^y/a^Zk - Nlh'^v:^^^/a^AkV^^ + Nlh^uvi^/a^AkVi, 
where 

Nlvh^Vi^kAkVi = 2Nlvhv^^kVH < N*u^h^Wi + akRl 
and, according to (16.121) and the inequality hlv^^l < NWi , 



-N*h^^kVy,AkV^, < N*Wi + il/3)akRl 



Therefore, 

/321 < N^h^^MZk + N*{1 + u^h^)Wi + {l/3)akRl + akR^ 
We can now specify h*: we take 

N*h* < 1 and h* < ho/2. 

Then, for h<h*, 

h2 < hy/TkZk + N*u^Wi + {l/3)fxakRt + (l/3)afci?2 + akR^, 



14/3 1 < h^Zk + N*u^Wi + N{v. 



70^ 
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+ {l/^)akRl + (l/3)/iafci?^ + N*atRl. (6-14) 
Estimating J4. By using Lemma [6.21 we easily see that 

I4/4I < N{v-^f + N*W^ + uWl'\N^k + N*h) Wil 

i 

where 

NvWl'^^kY, \vk^\ < NuWi + {l/?,)vakRl, 

i 

k,i 

It follows that 

I4/4I < N{v-f + N*uWi + {l/?>)uauRl. 

By combining this with (16.111) . (16.131) . and (16.141) . recalling that Zk > 
0, and coming back to (16.101) we conclude 

4J > (2afc - N^h^k)Zk 
- N{v-y - N*v^W^ + (z//3 - N;)akRl. (6.15) 
Now we specify = i^* by setting 

V* = 1 + ?,n; 

and finish the argument as in [13]. Namely, if 2ak — Nihy/a^ > for 
a k, then we can drop the term on the right in (I6.15P corresponding to 
this k because Zk > 0. However, if 2ak — Nih^fok < 0, then ^/o^ < Nh 
and \ak — Nihy/ak\ < Nh"^, whereas 

h'^Zk = h'^{v:^^AkVjg + fxv'jAkVij - uViAkVki) 

< N max W2 + NvWl'^ max Wl'^ 
~ Q Q 

< maxima + Nu^Wi < N(v-? + N*Wi. 
Q ^° 

This and (I6.15P yield (16.91) and the lemma is proved. 

Proof of Theorem 16.11 Fix a constant u according to Lemma 16.51 

and first assume that (to? 2:0) G Q. Then 

7/1 max < V^I^^{to,Xo) < N*{max \vij \ +max|t;i|), 

Q 02Q,i,j d2Q,t 



which by (16.41) yields similar estimate for 

max If 

After that (16.21) is immediate (cf. the end of the proof of Theorem 14. ip . 
Therefore, in the rest of the proof we assume that 
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Similarly, if (16.51) is violated, there is nothing to prove. Hence, we may 
assume that fl6.5p holds. Finally, we may assume that h < h*, where 
h* is taken from Lemma 16.51 and further reduced it if needed so as to 
satisfy h*y/h' < e. Indeed, ii h > h*, then in Q" 

\vii\ < 2(h*)~^ maxlvA. 

After justifying these additional assumptions which allow us to use 
the assertions of Lemmas 16.21 and 16.51 as long as h < h*, we construct 
functions f(t), ak(t,x), bk(t,x), c(t,x), /(p, t, x) = f(t,x) as in Sec- 
tion m to get (14. 2p and (14. 3p satisfied. Then, since (to, xq) G Q2 and 
( 16.31) and ( 16.51) are vahd and hy/u < e, by Lemma ESI at (to,Xo) we 
obtain 

^7oM^ (e"™^f5^t; + ak^kV + hdkV - (c + fCm)v + ^f) < 0. 

The fact that f, a^, and bk are limits of some r°, a^, and 6^, allows 
us to assert that Lemma [6751 holds with f, a^, and bk in place of r", a^, 
and 6^, respectively. Therefore, 

-N{v-^f + P,,^,{e-"'-f6Jv-{c + fcJv + ^f)<N*W^. (6.16) 

Here P^gfj_^6[^v > as right after (14.61) . Furthermore, 

Pioi^Aic + fc„,)v) = -(c + fCm)VJo^^•y + h, 

where Ii is a linear combination of products of two types: 

(i) f ~ or v~j times a difference operator applied to c times either w or a 
first-order difference operator applied to v- the second and third factors 
may be taken at a point different from {tQ,Xo), but their coefficients in 
the linear combination are dominated by a constant A^; 

(ii) Vi times a difference operator applied to c times either v or a 
first-order difference operator applied to v- these terms may be taken 
at a point different from {to,xo), the coefficients of these terms are 
dominated by a constant A^* (recall that u is entering P^p,p). 

Owing to Lemma 16.21 the absolute value of the linear combination 
of the products of type (i) is less than 

ATt;- i^3(max \v\ + Wl'^) < (t;" )' + Ar*(max \v\^ + Wi). 

The absolute value of the the linear combination of the products of 
type (ii) is clearly less than 

N*wl^'^(max \v\ + W^^^) < A^*(max \v\^ + Wi). 
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Now from the above estimates, f l6.16p . and the fact that V^gf^u > 
(f ~ )^ and c + fCm > A we conclude 

{X - N){v-f + P,,,,{^f) < N*im^x\v\' + W^). (6.17) 
Finally, obviously 

Pjof^A^f) > -iV*e(max \v,,\ + Umax \v,\) > -N*^ - ' 

and we infer from (I6.17P that 

{X-Ni){v-;)^ < N*{ma.x\v\'' + Wi)+N*i\ 

Q 

We set the constant in the statement of the theorem to be A'^i + 1 
and use Lemma W2\ to conclude that in Q° fl Q\o for any i,j 

\6,6M = \v,,\ < iVt;- (to, xo) < iV*(max \v\ + Wl^^) + N*^{to). 

Q 

This implies (16.21) in Q° n Q\o. On the remaining part of Qlo estimate 
(16.21) is obvious and the theorem is proved. 

7. Proof of Theorem 12.121 

We start with three auxiliary results. Everywhere in this section the 
assumptions of Theorem 12.121 are supposed to be satisfied. Recall that 
the set £ is introduced in ( ]2.20p . 

Lemma 7.1. For any function and li^h € Aq we have 

\KiM2<P{^)\ <4max(|Av,0(x)| : \k\ <do,xe (Aq + £) U {0}) 

+ 4max(|Aft,_40(x)| : do < \k\ < di,x e AqU {0}). (7.1) 
Proof. Obviously we may assume that h = 1. Next, observe that 
5i,iAhm = (l/2)[Ai,,,0(/i) + Ai,,,0(/2)] - (l/2)Ai,,,_,,0(O), 
5i,h'^i,ii0(O) = A^,/,0(/i), 5i,;,(5i,_z,0(O) = -A^,;,0(O). 
It follows that if li,l2 & C, then 

\Sh,iA,i2(P{^)\ < max{\Ah,eJ{x)\ : \k\ < do,x eCU{0}) 
+ max(|A,,,40(O)| : do < \k\ < di). 

We substitute here </>(?/ + ■) in place of and use that Ao + C D C since 
di > 2 and C = —C. Then we see that, if y G Aq U {0} and li, I2 G £, 
then 

\5i,iA,iJiy)\ < raax{\Ah,eJ{x)\ : \k\ <do,xe (Aq + £) U {0}) 

+ max(|A;,,40(x)| : do < \k\ < di,x E AqU {0}). (7.2) 
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In case = ^ + (2, k = Vi + V2 with d, (2, Vi, ^2 e £ and Ci 7^ C2, 
Ci 7^ -C2, Vi ^ V2, Vi -V2 either (i 7^ Vi and Ci ^ -Vi or Ci 7^ 
772 and Ci 7^ — '72- The second possibihty reduces to the first one by 
interchanging r]i and t]2. If the first possibility realizes, then we use the 
formula ^1,,,+,^ = Ti + to obtain 

'5i,Zi'^i,;20(O) = (^i,Ci'^i,C2 + ^i,Ci)(^i,'?i'^i,'?2 +Si,m)'Pi'^) 

= ^^l,C2^^1,r?2 0(Cl + ^71) + ^l,C2^1,m</'(Cl) + ^l,Cl^l,»)20(^l) + ^l,Cl^l,'?l</'(0)- 

Here Ci + ^71 ^ ^o^Ci)'?! ^ ^o? and G Aq U {0}. Therefore, we get 
(TO) from (17:21) . 

The remaining case that li E C and I2 = Vi + V2 with r/j as above is 
taken care of by setting (^1 = in the above calculations. The lemma 
is proved. 

Before stating the next lemma we remind the reader that the index 
k takes values in {±1, ±di}. 

Lemma 7.2. For any values of the arguments and s > we have 
2(1 

SQk, spk, stp) +Ko^qj: 

k k 

+Ks{J2\pk\ + \i^\ + r + s-')]. 

k 

Indeed, the expression in the brackets obviously is bigger than 

sup 

which in turn is bigger than 6q^ for each particular n = ±1, ±di. 
Below we use the notation T{e) and from Section O 

Lemma 7.3. Let 9 E (0, e^^) and e G (0, 1] be such that 

2diKoK/6 <l/2, K:=e{e + e)/{l-ee). (7.3) 
Let w, ip, pk be functions on Q and assume that 

m^x i^P^w)- < Omax'^l^Ah/M- (7-4) 

k 

Then in Q we have 

k 

where 

I = max[^Fi(f),Ah,e,w,pk,^) + Ks^{ \Pk\ + \^\ + (j)' + l)]- 

\k\<di 
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Proof. Set = J^ki^^h/k'^)'^ and observe that due to f l7.4p 

- e^^^^ > -9 max - max 
Q Q 

in Q. Hence {9 + e) maxQ > {e^^ — 9) maxg that is 

max$^ < Kmax^"*". (7.5) 

Q Q 

By (17. 5p and Lemma \T72\ with s'^cj), s~^Ah/f.w, s~^pk, and s~^ip in place 
of 0, qk,Pk, and ■?/', respectively, and = ^ we find that in Q 

< ^ [/ + KoK max $+1 . 

6 ^ Q 

Upon taking the maximums over Q of both parts and taking into 
account (17. 3p we get that < {Mi/ 6)1 in Q, which along with (17.50 
yield the result. The lemma is proved. 

Proof of Theorem I2.12[ Here k,i,j run through ±1, ...,±di. It 

is easy to see that one can find e = e{6,di,KQ) G (0,1] and = 
fi{6, di, Kq) > in such a way that the conditions: 

lQd\ii < lQdl{'iiifl'^ =:9< e'^ 

and (17.31) are satisfied. We choose and fix appropriate e and ii. 
If 

max {CPju)~ < 6lmaxy^ \^Ah/,^u\, (7.6) 
r(£)xQ Q ^—^ 

k 

then, by taking into account that u satisfies (12.11) in Q, from Lemma 
17.31 we obtain that \^Ah/f.u\ are bounded in Q by the right-hand side 
of (I2.2ip . By combining this with Lemma 17.11 we conclude that (12.2 ip 
is true in Qlo n Q2 (notice that (Aq + C) U {0} C Aq + Aq). Of course, 
(IMID is obvious on Qlo n d2Q. 

If (16. ip holds, then we get (I2.2ip from Theorem 16. 1[ In the remaining 
case both (17. 6p and (16.10 are violated and 



^i^ :=^max^|^A 



2\l/2 



< max (^P^m) < (3/i)^/^ max ( ^ \^6h,ej6h/M'^) 
< (3/i)^/^max \^5h,e.5h,eM + {^f^Y^^ raaxS^ \^5h,e,5h,eM- 



''2 

72; 



In light of Lemma 17.11 the last maximum over Q2 is less than d>d\H. 
Hence 

9H <N max le^h^/h/.w! + (l/2)^i/, H < N max l^Sh/Jh/M 
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and we can finish the proof (12.211) as few times before. The theorem is 
proved. 

8. Proof of Theorem 12.141 

In the following lemma the assumption that sup„ > 5 is not used. 
All other assumptions of Theorem 12.141 are supposed to hold. We use 
notation f l2.10p and the notation from Section [6] with + 1 in place of 

di and 5i = 6h,A, = Ah,,e^. Here we take 

/i = 

and show how to choose u = z/(A, e, di, K^) in Lemma 18. 1[ 

Lemma 8.1. In Q\q for 7 G r(£:) we have 

(V7.A;,^,,^m)- < iV*e™^(^+^)(l + max|e(_)5,^,,^n| 
^ Qd 

+ max + max(^(_)Aft^,^^M)"), (8.1) 

where N* = N*{\,ho,e,di,Ks). 

Proof. As many times before, if {to,XQ) G diQ, there is nothing to 
prove. Therefore, we assume that ito,Xo) G Q?. We may also assume 
that at {to, xq) 

i i 

Then the operator P-yoOi^ respects the maximum principle for hv < Is 
(see the proof of Lemma 16.21 and recall that rj < h). 
Then as in the proof of Lemma 16.51 we obtain 

4J = AP^^ouiak^kV + bkdkv) > au^Rl + 4J4 (8.2) 

if h < /?-o/2, where this time 

h = 2v~^-foi{6ibk)vki + v~^-foi{Aibk)vk - i^Vi{Sibk)Th^/,Vk. 

Since 

v-,loi\{Sibk)vki\ < N*v-^^Y1 < (A/8)(t;-)' + N*akJ2^li 

i i 

we get from (18. 2p that, for u = z/*(A, e, di, K^), hu* < and h < hojl, 

J>-[\/2){v-;)'-N*W,. (8.3) 

The rest is just a repetition of a part of the proof of Theorem 16.11 
with obvious and big simplifications. The lemma is proved. 

There is almost nothing else to do to finish the proof of Theorem l2.14[ 
Indeed, (18.11) with di in place oi di + 1 yields the first estimate in (12.221) 
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as in the proof of Theorem 12.121 After getting estimates for lA/i^^^-ul 
the estimate of (A^ ^m)^ follows immediately from fl8.ll) . The theorem 
is proved. 

9. Comments on the operators having form (12. 7p 

We know (see, for instance, [^) that if an operator L having form 
(12.71) admits an approximation with operators Sh of the form (12.41) 
respecting the maximum principle with hB in place of B and Span B = 
M.'^, then necessarily 

Lu = O'k^l.il.u^ixi + bkil-U^i 

with some a^, bk > and ik E BU {—B). A way to find such represen- 
tations for (i = 2 and given a^^ is suggested in [Ij. 

Next natural issue is related to the smoothness of if we are 

given that a*-' are smooth. Recall that in Assumption 12.21 we need 
to be at least Lipschitz continuous. Of course, this problem disappears 
if a*-' are constant. 

It is an easy and probably well-known fact that if (a*-') = (a*-' (t, x)) 
is uniformly bounded and uniformly elliptic, then one can find di and 
Aq for which can be chosen strictly positive and as smooth as a*-^ are. 
The proof of this can be obtained from the fact that if we are given 
a closed convex polyhedron then every point in the relative interior 
can be written as a convex combination of the extreme points with 
the coefficients > which are infinitely differentiable functions of the 
point. By replacing L with L + e^A one can approximate a possibly 
degenerate operator L with uniformly nondegenerate ones, so that there 
always exist a sequence of operators of the form (12. 7p approximating L. 
Notice, however, that generally the set Aq changes with e. Nevertheless, 
one knows how to estimate the difference of solutions corresponding to 
L + e^A and L (see, for instance [9]) and between the solutions of 
the corresponding finite-difference approximations (see, for instance. 
Theorem 5.6 of [13j or Remark 12. 6p . 

Another generic example is given by the so-called diagonally dom- 
inant matrices. For instance, take d = 2 and assume that 6 = 0, a*-' 
are twice continuously differentiable with respect to x, a^^ = a^^, and 
|a^^| < s, where s = a^^ A a^^. 

Set K = 1/3 and take an infinitely differentiable, even, and convex 
function ilj{t) on M such that ipiu) = lul fo^' lul ^ Introduce 

g = a^^s-\ h = silj{g), 2d^'^'^ = h ± 2d'' = a'' - h, 
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where i = 1,2 and 0-0 ^ := 0. For other values of i,j = ±1, ±2 define 
a^^ so that 

and set 

£j = e[j|sign j, i^j = £i + ij^ 
where ei, 62 are the basis vectors. Then simple manipulations yield 

ALu = a'^iijiljU^k^r. 

We now show that not only L admits a representation as the sum of 
second order directional derivatives with the directions independent of 
t, X but also Vd^ are Lipschitz continuous in x. By the way, observe 
that obviously a*-' > 0. 

We are going to use that nonnegativc and twice continuously diffcr- 
entiable functions are the squares of Lipschitz continuous functions. In 
particular, a" and, consequently, a^^ A a^^ are the squares of Lipschitz 
continuous functions, |a"| < N\/a^ and \sx\ < N^^. Furthermore, 
a^^ ± a^^ is nonnegative and twice continuously diffcrcntiablc. Hence, 
it is the square of a Lipschitz continuous function. In particular, 

lal"- ± af\ < NVa^^±d^, \al^\ < N^/^, \al^\ < N^/^. 

and recalling that |a^^| < s we find 

\gx\ < K^ls'^ + \a^^\ ■ \sx\s~^ < Ns'^/^, \h^\ < Ny/^. 

Next, the function := ip{y) + y is smooth and nonnegative. 

Therefore 

Similar estimate holds for |a^'~^| and |a^ | if i 7^ j. 

On the set where |a^^| > ks, we have 2a}^ — a^^ — \a^^\, so that by 
the above 

Finally, on the set where |a^^| < 2ks it holds that h < 2ks, a^^ — h > 
na^^, and {a,]}] < N\/a^ < N\/dM. Similarly we get what we need 
for a^^ and the remaining a". 
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